Rectification and Flux Reversals for Vortices Interacting with Triangular Traps 
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We simulate vortices in superconductors interacting with two-dimensional arrays of triangular traps. 
We find that, upon application of an ac drive, a net dc flow can occur which shows current reversals 
with increasing ac drive amplitude for certain vortex densities, in agreement with recent experiments 
and theoretical predictions. We identify the vortex dynamics responsible for the different rectification 
regimes. We also predict the occurrence of a novel transverse rectification effect in which a dc flow 
appears that is transverse to the direction of the applied ac drive. 
PACS numbers; 74.25. Qt 



Vortices interacting with periodic pinning structures 
have been attracting considerable attention since they are 
an ideal system in which to study the static and dynami- 
cal behaviors of collectively interacting particles coupled 
to periodic substrates. In these systems, a variety of com- 
mensurability effects and novel vortex crystals can occur 
as a function of the ratio of the number of vortices to the 
number of pinning sites [1~3]. The vortex crystals are 
composed of one or more vortices trapped at each pin 
with any remaining vortices located in the interstitial re- 
gions between the pinning sites. A variety of dynamical 
flow phases occur when a drive is applied, including the 
flow of interstitial vortices between the pinning sites [4,5] 
as well as channeling along rows of pinning [4] . Much of 
the physics of vortices interacting with periodic pinning 
is also observed for repulsively interacting colloids in pe- 
riodic optical trap arrays [6-8]. In addition to the basic 
science issues, vortices interacting with periodic pinning 
arrays are relevant to possible technological applications 
of superconductors, such as critical current enhancement 
or controlled motion of flux for new types of devices. 

Several methods have been proposed for using periodic 
pinning or controlled disorder in superconductors to cre- 
ate vortex ratchets, rectifiers, and logic devices [9-13]. 
As in general ratchet systems, which can be determin- 
istic or stochastic in nature [14], a vortex ratchet trans- 
forms an ac input into a dc response. Vortex ratchets can 
be created via asymmetric periodic pinning lines [9,10], 
asymmetric channels, or the asymmetry introduced by 
multiple ac drives in a system with a symmetric sub- 
strate [11,12]. Ratchets constructed of periodic two- 
dimensional (2D) arrays of asymmetric pinning sites have 
also been proposed [13]. Recently, both positive and neg- 
ative vortex rectification have been experimentally real- 
ized in periodic arrays of triangular pinning sites [15]. In 
Ref. [15], the triangles are arranged in a square lattice 
with the tips of the triangles oriented in the +y direc- 
tion. When an ac force is applied in the y direction at low 
matching fields, rectification of the vortex motion in the 
+y direction occurs over a range of ac amplitudes, with a 
peak rectification at a particular amplitude. For higher 
matching fields, there is a —y rectification at lower ac 
drives, followed by a +y rectification at higher ac drives. 



This change in the rectification direction is explained in 
Ref. [15] as originating from the separate motion of inter- 
stitial vortices at low drives, giving —y rectification, and 
the motion of the pinned vortices at higher drives, pro- 
ducing +y rectification. It is not clear how the motions 
of these two vortex species can be fully separated, espe- 
cially in the presence of thermal fluctuations, and thus 
a clearer picture of the vortex dynamics producing the 
rectification is needed. 

In this work we present simulations of vortices inter- 
acting with 2D square arrays of triangular pinning sites 
for parameters corresponding to the recent experiments 
of Ref. [15]. We find a +y rectification at low ratios n 
of the number of vortices to the number of pinning sites, 
with a maximum +y dc flux as a function of ac amplitude. 
When interstitial vortices are present, the initial rectifica- 
tion is in the —y direction, and is followed by a crossover 
to +y rectification at higher ac amplitude, as seen in ex- 
periments. We find that for large ac amplitudes, when 
the motion of the vortex lattice becomes elastic, the rec- 
tification switches to the —y direction, which is explained 
in terms of channeling effects. We also predict a novel 
transverse ratchet effect where a dc motion is generated 
in the direction transverse to the applied ac drive. All of 
the effects we describe here should also occur for repul- 
sively interacting colloidal particles interacting with 2D 
triangular traps. 

We consider a thin film superconductor containing a 
square array of Np = 64 triangular pinning sites with 
periodic boundary conditions in the x and y-directions. 
We add Ny = nNp vortices to the sample, each of which 
obeys the overdamped equation of motion 



dt 



f„ + f^c + fT 



(1) 



Here the vortex-vortex interaction force, appropri- 
ate for a thin film superconductor, is ~ 

J2j^ii^o/ f^o^^)"^ /'''ij : where $o is the elementary flux 
quantum and A is the thin film screening length [16]. We 
use a fast summation method to evaluate this long-range 
interaction [17] . The pinning force fp comes from a square 
array of equilateral triangles with one vertex pointing in 
the +y direction, as shown in Fig. 1. Every triangular 
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FIG. 1. Pinning sites (triangles) and vortex positions (dark 
circles) at the matching fields n = (a) 1, (b) 2, (c) 3, and (d) 
4. 



pin is modeled as three half-parabolas of equal strength, 
each of which attracts the vortex to a line passing through 
the center of the triangle and parallel to one of the sides. 
The pinning force is cut off at the edge of the triangle. 
The applied ac force {ac = Asin{ujt)r acts in the di- 
rection perpendicular to an applied ac current. Here we 
consider either J^q = Asm{ujt)y or f%(y = Asm{uit)x, 
with no mixtures of ac drives. Temperature is modeled 
as random thermal kicks with the property < fr(i) >= 
and < tT{t)tT{t') >= 2r]kBTS{t - t'). 

We match our parameters to those of the experiment 
in Ref. [15], performed in Nb films with 77 = 1.4 x 10^^^ 
Ns/m. We take T/T^ = 0.98, giving /t = 0.46/o, where 
/o = 0g/27r/ioA3 = 1.09 x IQ-^ N/m. At this tem- 
perature, the London penetration depth A = A(0)/(1 — 
(T/rc)2)i/2 = 368 nm, so our pin spacing is 2.09A in the 
X direction and 2.03A in the y direction, and the pins 
are 1.68A on a side. In the experiment, each pin held a 
maximum of three vortices, so to match this we set the 
pinning strength to fp — 1.05/o. We fix = 78 kHz, 
higher than experiment due to the limitation of simula- 
tion time; however, since the behavior of the system is 
controlled by the ratio A/uj, we compensate by taking 
A larger than in experiment. We note that we find the 
same generic behaviors for other parameters, including a 
pinning array that is triangular rather than square. 

In Fig. 1 we show the vortex positions and pin locations 
at yl = for a system with a square array of triangular 
pins at fillings n = 1, 2, 3, and 4. A global symmetry 
breaking occurs aX n — 2 between the two possible ar- 
rangements of the vortices allowed by the square pinning 
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FIG. 2. Net dc velocity <Vy >in units of m/s as a function 
of applied ac current A in units of /o for n = (a) 1, (b) 2, (c) 
3, (d) 4, (e) 5, and (f) 6. Inset to (a): dc depinning force \fdp\ 
for -f y (circles) and — y (x's) directions. 

lattice. All of the vortices sit in the pins for n < 3, but 
for n > A the vortex-vortex interaction is strong enough 
that some vortices move to the interstitial regions. 

We first consider the dc depinning force fdp in the pos- 
itive and negative y directions by applying a dc driving 
force of increasing magnitude and measuring the average 
vortex velocity <Vy>. In the inset to Fig. 2(a) we plot 
\fdp\ as a function of n. At n = 1 the depinning is sym- 
metric in the ±j/ directions, as expected given our model 
for the pinning. For n — 2, l/^"] is twice as large as 

l/dffl' "^^^ ^^"^^ ih.aX the vortices align vertically for 

+y driving, and horizontally for —y driving. For n — Z 
and 4, \f2p\ is 25% larger than |/^^|. Here, the vortex 
at the top of each pin acts to assist in the depinning of 
the vortices along the bottom of the pin for — y driving. 
For rt > 4, the depinning is dominated by the intersti- 
tial vortices and the asymmetry in the depinning is lost. 
The depinning asymmetry at n = 2 to 4 should be large 
enough to be observable experimentally. Depinning in 
the ±x directions is symmetric for all n. Asymmetric 
depinning at noninteger n was observed in Rcf. [18]. 

We next consider the effect of an applied ac drive f^^j . 
We monitor the net dc velocity < Vy > aX fixed lo and 
sweep the ac amplitude A, averaging <Vy > for 20 peri- 
ods at each increment of A. In Fig. 2(a) we plot <Vy > 
vs A for n — 1. For A < 0.9, there is no net flow in x or y, 
indicating that the vortices are still pinned. For A = 0.9 
a finite +y dc velocity appears with no net x dc velocity, 
corresponding to a +y rectification. As A increases fur- 
ther, < Vy > increases, reaching a maximum at A = 0.99 
and then dropping back to zero at ^ = 1.05. This behav- 
ior matches that seen in the experiments of Ref. [15]. We 
find no steps in the +y rectification regime, in agreement 
with the experiments, but in disagreement with recent 
T — simulations of a very small system for asymmetric 
pinning sites that predicted steps in the velocity vs 
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FIG. 3. Pins (triangles), vortices (dots), and vortex trajec- 
tories during 1/4 of the driving period (lines) for: (a) N = 2, 
A = 1.0, —y portion of drive cycle, (b) +y portion of same 
cycle; (c) A'^ = 4, j4 = 0.35, —y portion of drive cycle, (d) +y 
portion of same cycle; (e) A'^ = 4, A = 0.75, —y portion of 
drive cycle, (f) +y portion of same cycle. 

/ac curve [19]. The transport does not occur in an orga- 
nized step like flow, but is instead a stochastic process, 
as we will show in more detail in Fig. 3. In Fig. 2(b) and 
Fig. 2(c) we plot < > for n = 2 and 3, respectively. 
These curves are similar to the n = 1 case, with an ini- 
tial pinned phase followed by +y rectification that goes 
through a maximum with increasing A. 

In Fig. 2(d) we show < Vy > for the case of n = 4 when 
interstitial vortices are present, as seen in Fig. 1(d). Here 
there is an initial negative rectification regime in the —y 
direction, followed by a crossover to a +y rectification 
as seen in experiments [15]. In Fig. 2(e,f) we plot the 
cases for n = 5 and n — 6, respectively, which show a 
similar rectification to that in Fig. 2(d). In all of these 
combination of interstitial and pinned vortices 
are present. As n is increased further, the ratcheting ef- 
fects are gradually reduced and become indistinguishable 



from thermal noise. 

For n = 1 in Fig. 2(a), at higher A > 1.05, < Vy > be- 
comes negative, indicating —y rectification at high drives. 
This reversal occurs for all fillings but is most pronounced 
for n = 1 and n = 4, when the onset of the reversal falls 
at lower values of A. We predict that this effect should 
be easily observable in experiment, although it was not 
reported in Ref. [15]. The high-drive —y rectification 
originates when the — y vortex motion occurs through in- 
terstitial channels while the +y vortex motion becomes 
localized on top of a column of pinning, giving a larger 
drag during the +y portion of the cycle. 

To understand the role of interstitial and pinned vor- 
tices in the rectification process, we examine the dynam- 
ics of the vortices. In Fig. 3(a,b) we illustrate a typical 
example of vortex motion leading to +y rectification at 
n = 2 and A = 1.0. Fig. 3(a) shows the vortex trajecto- 
ries during 1 /4 of the drive cycle when the drive is in the 
—y direction. At these temperatures, none of the vortices 
remain pinned. The vortex motion is uncoordinated and 
the vortices wander stochastically as they pass through 
the pinning. In contrast, during the +y portion of the 
cycle shown in Fig. 3(b), the vortices are channeled both 
within and between the columns of pinning sites, and 
there is little wandering in the x direction. This more 
focused motion leads to a greater overall displacement 
in the +y direction, and a net +y rectification. Similar 
motion occurs for n = 1 and n — 3, except that for n = 1 
all of the motion is confined to the pinning channel and 
there is no interstitial motion, so the net rectification is 
correspondingly reduced as seen in Fig. 2(a). 

At n = 4 and above, when interstitial vortices are 
present, —y rectification occurs at low drives, as shown 
in Fig. 2(d-f). In Fig. 3(c-d) we illustrate an example 
of vortex motion in the —y rectification regime for n = 4 
and A — 0.35. Fig. 3(c) shows that the vortices meander 
moderately in the x direction during the —y portion of 
the drive cycle. Note that none of the vortices remain 
pinned; instead, all of the vortices are participating in 
the motion, with pinned and interstitial vortices switch- 
ing places frequently. During the +y portion of the drive 
cycle shown in Fig. 3(d), the drive is too low to overcome 
the pinning force and allow the formation of channels of 
the type seen in Fig. 3(b). Instead, the +y vortex fiow 
is strongly diverted into the x direction by the pins, re- 
ducing the net motion in the +y direction, and resulting 
in an overall —y rectification. As A is further increased, 
the vortices begin to overcome the pinning in the +y di- 
rection and form channels during the +y portion of the 
cycle, as illustrated in Fig. 3(f) for rt = 4 and A = 0.75. 
This gives +y rectification just as in Fig. 3(a-b). As n in- 
creases above 4, the channels become more closely packed 
and the amount of vortex wandering in the x direction 
during the —y portion of the cycle decreases, progres- 
sively lowering the net +y rectification. The onset of the 
—y and +y rectification phases also drops to 
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FIG. 4. Transverse rectification for /ac- (*) < Vx > for 

n = 2; (b) corresponding < Vy > with —y rectification, (c) 
< 14 > for n = 4; (b) corresponding < Vy > with +3/ rectifi- 
cation. 

lower A with increasing n since the effectiveness of the 
pinning decreases with increasing vortex density. 

In Fig. 4 we iUustrate < Vx > and < Vy > when an 
ac drive f%Q is apphed in the x direction for n = 2 and 
4. The net x velocity is zero but we find a transverse 
rectification in the y direction: —y for n = 2 and +y for 
n = 4. The magnitude of the rectification is comparable 
to that seen for ,fj^(^, and it should thus be experimentally 
observable. The transverse ratchet efi^ect is produced by 
the interaction of the vortices with the pinning sites. For 
n = 2, the vortices channel along the rows of pins and 
are defiected downward by the pinning sites, producing 
—y rectification. Similar motion occurs for n = 3. For 
n = 4 to 6; the interstitial vortices dominate the motion 
through channels that pass along the tips of the triangles, 
although all vortices are being pinned and depinned. The 
interstitial vortices are deflected upward by the vortices 
inside the pins, producing +y rectification. Since the 
transverse motion does not require breaking of the x di- 
rection symmetry of the system, it can also be observed 
for a dc a; direction drive. 

In conclusion, we have conducted simulations of vor- 
tices in thin film superconductors interacting with trian- 
gular pinning sites for parameters relevant to recent ex- 
periments. We show that both a +y and —y rectification 
can occur at depinning depending on whether interstitial 
vortices are present, in agreement with experiment. We 
find dc depinning anisotropy for n = 2 to 4, but observe 
that even when the dc depinning is not anisotropic, rec- 
tification can still occur due to the fact that the vortex 
dynamics differs under ac and dc drives. Rectification in 
the +y direction occurs when the ac drive overcomes the 
pinning strength and vortex channels form that flow in 
the +y direction. When interstitial vortices are present 
at zero ac drive, an initial —y rectiflcation occurs at lower 
ac drives due to scattering of the +y vortex motion by the 
pinning sites. This occurs for n = 4 and above, in agree- 



ment with experiments. We also predict that a novel 
transverse ratchet effect should occur when the ac drive 
is applied in the x direction. Here, the rectification is 
in the y direction. Our results should be generic to any 
type of repulsively interacting particles moving through 
triangular traps and may have potential applications for 
species fractionalization in colloidal systems. 
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